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In order to test different models proposed for some states discovered recently in the char- 
monium mass range that do not fit into the pattern predicted by the conventional quark odel, 
we derive predictions for the spectrum within the hadro-charmonium picture, the tetraquark 
picture as well as the hadronic molecular approach. We exploit heavy quark spin symmetry 
for the hadro-charmonium and hadronic molecule scenarios. The patterns that emerge from 
the different models turn out to be quite distinct. For example, only within the hadro- 
charmonium picture a pseudoscalar state emerges that is lighter than the Y (4260). Possible 
discovery channels of these additional states are discussed. 


PACS numbers: 14.40.Rt, 13.75.Lb, 13.20.Gd 


I. INTRODUCTION 


In the past decade a number of states in the charmonium and bottomonium sectors have been 
discovered with properties in conflict with the quark models for mesons which were quite successful 
in describing the low-lying heavy quarkonium states as QQJ^ = c,b) states, such as the classical 
Godfrey-Isgur quark model [l| and Gornell potential model 2-^. The newly observed structures in 


the heavy quarkonium mass region include, am ong many others, the X(3872) E(4260) B and 


the charged states Zc(3900) 
and the Zfe(10610), Zfe(10650) 


, Zc(4020) 1 Ch 13| and Zc(4430) 14l4l6| in the charmonium sector 


171 ] in the bottomonium sector (see, e.g. the mini review on heavy 


‘Electronic address: eleven®ihep.ac.cn 
^Electronic address: fkguo@hiskp.uni-bonn.de 
^Electronic address: c.hanhart@fz-juelich.de 
^Electronic address: q.wang@fz-juelich.de 
^Electronic address: zhaoq@ihep.ac.cn 













2 


quarkonium spectroscopy in the Review of Particle Physics by the PDG [3]). Because of the failure 
in describing these structures simply as QQ mesons, they are considered as candidates of exotic 
hadrons^. Various models were proposed to explain these structures, see e.g. Ref. H . However, 
no consensus for almost any of these states has been achieved. It is thus of utmost importance 
to scrutinize these models, and manifest their distinct predictions. Because the near-threshold 
narrow structures in the continuum channel of the open-flavor meson pair cannot be explained by 
just a threshold cusp [21[, we will assume that all the states to be discussed correspond to physical 


states. The dynamical structure of these states are being investigated in several different models. 
The purpose of this work is to present distinct predictions of several commonly discussed models. 
The scenarios to be considered include the tetraquark model, the hadro-charmonium model, and 
hadronic molecules. A brief description of each model including a list of relevant references and a 
short description of the underlying assumptions will be presented in the corresponding sections. 

One has to keep in mind that different components mix with each other whenever such a mixing 
is not forbidden. Thus, when we present predictions for a certain scenario for a given state, we 
always mean that we are looking for the consequences that arise if that scenario is assumed to be 
dominant. 

The key feature we are going to exploit in this work is heavy quark spin symmetry (HQSS). 
It is an approximate symmetry and becomes exact in the limit of infinitely heavy quarks. It 
arises because the spin-dependent quark-gluon coupling in quantum chromodynamics (QCD) is 
proportional to the magnetic moment of the heavy quark, and vanishes in the heavy quark limit. 
As a result, it leads to spin multiplets. Within a multiplet, the masses are degenerate in the 
heavy quark limit, and the mass splittings depend on the dynamics of the model for finite quark 
masses. As will be discussed, the multiplet structure differs in different scenarios, and thus provides 
invaluable information. We will compare the spectroscopy predicted by each of these three scenarios 
in the following sections. We close with a summary. 

II. HADRO-QUARKONIUM 
A. Assumptions 


A hadro-quarkonium is a system with a compact heavy quarkonium embedded inside a cloud of 
light hadronic matter [2, 23]. This scenario was motivated by the fact that several charmonium- 


^ In fact, in the seminal quark model paper [l^, Gell-Mann already mentioned the possibility of multiquark states. 
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like states were only observed in final states of a specific charmonium and lig 
are the y(4260) discovered in J/ipTr-K Q], the Zc(4430) discovered in -^'tt 


y(4660) observed in ip'mr 


M 


it hadrons. Examples 


24l |. the y(4360) and 


The recent BESIII observation of similar cross sections for 


J/7r~ and /icvr"*‘ 7 r“ at 4.26 GeV and 4.36 GeV in e'’'e“ collisions [?, 11| stimulated Li and 
Voloshin to generalize the hadro-charmonium model to include HQSS breaking and describe the 
V(4260) and y(4360) as a mixture of two hadro-charmonia [^ : 

Y {4:260) = cos 9 — sin 6 ipi, Y {4360) = sin 9 + cos 6 ^jJl , (1) 

where 'tpi ~ ( 1 "' )ccG)(0 ^)qq and 'i/’s ~ (1 )cc® (0''"'’)gg are the wave functions of the = 1 
hadro-charmonia with a 1"' and 1 cc core charmonium, respectively. It was argued in Ref. 


27| 


that t /’3 contains predominantly a ip' (rather than the ground state J/ip) and that V’l contains 
predominately an he- The decays into the J/ipTT~^Tr~ then occur through de-exciting pj' to J/ip in 
the light hadronic matter. The strength is controlled by the so-called chromo-polarizibility ol:^' jj^^, 
(see Ref. Q)- 

If we assume that the Jl'(p'K^'n~ events seen in e+e“ collisions at energies around 4.26 and 
4.36 GeV are mainly from the decays of V(4260) and V(4360), the BESIII data imply that 


r(y(4260) —)■ J/'ip'K'n) is similar to 


angle as large as 0 ~ 40° is needed 


'(y(4360) —)• J/'ip'KTr). To achieve this pattern a mixing 


27j. Such a large angle translates into a small mass difference 


between the pJi and ^3 hadro-charmonia since 


tan( 20 ) = 


2m 


13 


TTlip^ 771-03 


( 2 ) 


which can be obtained from 


I "^Y(4260) 0 

0 l«y(4360) 


I 


cos y — sm 
sin 9 cos 9 



cos 9 sin 9 
— sin 9 cos 9 


( 3 ) 


77Z01 771]^3 

777i3 77703 

A mixing angle of around 40° leads to « 4.30 GeV, « 4.32 GeV and a mixing amplitude 
mi 3 50 MeV. 


B. Consequences 


In addition to the interference patterns in the line shapes discussed in Ref. 


271], what else does 


the proposal of the Y (4260) and Y (4360) as mixed hadro-charmonia imply? As mentioned in the 
Introduction, HQSS is quite useful in this respect. The binding force between the charmonium core 
and the surrounding light hadronic matter is due to the exchange of soft gluons. Because both the 
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charmonium and light hadronic matter are color singlets, at least two gluons need to be exchanged. 
The leading order (LO) interaction is due to exchanging two chromo-electric gluons H . It is 
important to notice that the heavy (anti-)quark spin decouples from such an interaction. Therefore, 
the LO interaction between the light hadronic matter with a particular heavy quarkonium, [QQ], 
should be the same as that for the spin partner of that [QQ] state. This means that a hadro- 
quarkonium should have spin partner(s) just as its core heavy quarkonium does. 

Exchanging one chromo-electric and one chromo-magnetic gluon provides a P and CP odd force, 
and thus induces the mixing between two hadro-quarkonia with core heavy quarkonia of opposite 
P and CP. The mixing between the "03 and 'ipi states, which contain the 'll}' (P = — ,CP = -|-) 
and the he (P = -|-,CP = —), respectively, is such an example. Exchanging two chromo-magnetic 
gluons provides one source for splitting the masses for one spin multiplet of hadro-quarkonia. It 
is suppressed by ©(Ag^j^/mg) in comparison with the LO interaction, and gives a tiny correction 

4% for hadro-charmonium and ~ 1% for hadro-bottomonium) to the potential energy and thus 
to the mass of the hadro-quarkonium. The mass splitting between the spin partners within the 
same multiplet of hadro-quarkonia is therefore given approximately by the mass splitting between 
the core heavy quarkonia with the next-to-leading spin symmetry violation controlled by mixing 


analogous to the one discussed above for the y(4260) and the y(4360). In fact, in Ref. [30|, 


analogous to the question of interest, HQSS has been used to predict that the y(4660) as a 
"0^0 (980) bound state 3^ which may be regarded as a specific example of hadro-charmonium has 
a spin partner: an ry(,/o(980) bound state, with mass of around 4616 MeV. 

Erom the discussion above it follows that the "03 state has a spin partner characterized by the 
same light quark cloud with the xjj' in the core replaced by the ry).. We may call this state ryc[^c]- It 
has = 0 t and a mass of around 


m ^3 — {m^i — rrirj'J ~ 4.25 GeV. (4) 

Similarly, replacing the he in the V’l state by any of the Xcj states leads to three spin partners of 
"01. The quantum numbers of these states composed of (0 ^)qq are = J Their 

masses are listed in Table |T] and shown in Eig. [T]for illustration. We notice that there are two states 
with = 0 ^ in analogy to 1 . As in the vector channel the small mass difference of the pure 

spin states of about 40 MeV introduces a sizable mixing. The relation of the mixing amplitude in 
this case can be related to that in the vector case by constructing an CP-odd operator with HQSS 
breaking 

Omixing = \{x^ ■ ^J') + h-c- = xloh'c + , (5) 





5 


Composition 

Label 


Mass (GeV) 

■4' O {0++)qq 

V'3 

1 — 

4.30 

ic ® (0++),,- 

heWe] 

0-+ 

4.25 

he® (0-+),, 

jjjl 

1— 

4.32 

XcO ® (0 ^)qq 

»7c[Xco] 

0-+ 

4.21 

Xcl ® (0”+)qg 

Vcl [Xcl] 

1-+ 

4.31 

Xc2 ® (0"+)99 

??c2[Xc2] 

2-+ 

4.35 


TABLE I: Masses and quantum numbers of hadro-charmonia in the spin multiplets of the ipi ^ he 0 (0 
and ipz '4'' ® states. 

4.40 

4.35 

^ 4.30 
> 
u 

a 4.25 
I 4.20 
4.15 
4.10 




F(4360) 

7c2(4350) 


77X432U) ,7^,(4310) 


■ F(4260) 

. 

: 

77X4140) 

: 1““ 

0-+ 1-+ 2“+ ; 




FIG. 1: The spectroscopy of the spin partners of F(4260) and y(4360) based on the assumption that they 


are mixed states of two 1 hadro-charmonia as proposed in Ref. 


27|. The dotted lines in the vector and 


pseudoscalar sectors represent the masses of the unmixed states. The masses of y(4360) and E(4260) were 
used as input for the analysis and are shown as dashed lines. 


where a are the Pauli matrices and (...) takes trace in the spinor space. The fields J' and 


annihilates the S-wave and P-wave charmonium states, respectively, and are given explicitly by 




ml, 


j'= a+ X* = (-XS - + K • (6) 

Thus the mixing amplitude is up to corrections of order (Aqcd/^tIc)^, with mis ~ 50 MeV. 

We can calculate the eigenvalues of the mass matrix ^ 


"^ 7 ,^ 4320 ) ~ 4.32 GeV , and m^^( 4 i 4 o) « 4.14 GeV , 


(7) 


^ We thank M. Voloshin for pointing this out. 

® Because these states have pseudoscalar quantum numbers, we name them as ryefmass). 
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and the mixing angle is Oq ~ —38° with 6q defined via 


rye(4320) = cos 00 ryc[ry'] - sin 0o hc[Xco], rye(4140) = sin 0o hc[^c] + cos 0o ??c[Xco] 


( 8 ) 


Please observe that among the states predicted there is a state with exotic quantum numbers 
1 it appears to be a quite robust consequence of the hadro-charmonium scenario proposed in 
Ref. that an rjd state exists with a mass between those of y(4360) and y(4260). 

Since all the considered hadro-charmonium states are above the corresponding thresholds for 
the decay into the core charmonium and two pions, we expect that they decay easily through 
dissociating the light hadronic matter into two pions, in complete analogy to Y (4260) and Y (4360) 
that were observed in final states consisting of a charmonium and two pions. This mechanism will 
introduce a width of ~ 100 MeV, in the ball park of the widths of the Y (4260) and the Y (4360), for 
each of them. Considering that both the widths of the ry(, and the XcO are about 10 MeV 1^, much 


larger than those of the ijj' and the he, the predicted states 7yc(4140) and 7yc(4320) can also decay 
via the decays of the Ty(. and Yro- Hadro-quarkonia can also decay into open flavor heavy meson 


and antiheavy meson pairs 


34l |. but it is natural to expect the first mechanism to be dominant. 


In the above, we have argued that if the Y (4260) and the Y (4360) are mixed hadro-charmonium 
states, then it is very likely that they have spin partners as shown in Fig. [T] Thus searching for 
these partners can provide valuable information on the nature of the y(4260) and y(4360). In 
which processes should they be searched for? As for the two pseudoscalar states ?7c(4140) and 
ryc(4320), because of the probably sizeable mixing, both of them decay into the ryc^Trvr and XcO'^’T- 
Being pseudoscalars, they cannot be produced directly in e^e“ collisions. One way of searching 
for them is to measure the rye invariant mass distribution for the decays —>• K^ric'^TT^TT~ 

as suggested in Ref. H for searching for the spin partner of the y(4660). Another possible way 
of searching for the 7yc(4140) and Tyc(4320) is to study the radiative decays of the y(4260) and 
y(4360). Because the branching fraction for ijj' —?> yxco is 10% [l^ . two orders of magnitude larger 
than that for ijj' yry^,, both the 7yc(4140) and 7yc(4320) states can be produced through the decays 
of the ips components of the y(4260) and y(4360) into their ryc[Xco] components. Hence, one may 
search for these two states in the process e“*"e“ —>■ "yXcoV at the center-of-mass energies around the 
masses of the y(4260) and the y(4360). 

The exotic state ryci(4310) and the state Tyc2(4350) can be searched for in analogous processes in 
the decays of Y (4360) with the XcO in the final state replaced by the Xci and Xc 2 , respectively. This 
kind of measurements may be performed at BESIII and a future high-luminosity super tau-charm 
factory. 
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If there are hadro-charmonium states, it is natural that the analogous hadro-bottomonium 
states should exist as well 2^. The interaction strength between the heavy quarkonium and 


the light hadronic matter is dictated by the chromo-polarizibility which is a matrix element for 
the propagation of a color-octet QQ pair [29|, and thus depends on the wave functions of the 


heavy quarkonia involved. Numerically it was found that the off-diagonal chromo-polarizability 
for the transition T' —)• Ttttt is a factor of 3 smaller than that for the charmonium analogue 
28l|.^ Therefore, there is no flavor symmetry connecting the hadro-bottomnium 


ijj' J/lp'K'K 


to the hadro-charmonium. However, because the mixing is induced by exchanging one chromo¬ 
electric and one chromo-magnetic gluon, one naively expects that the mixing amplitude for hadro- 
bottomonium states is much smaller than that for the hadro-charmonium states, and roughly scales 
down by a factor of mc/ruf,. If the unmixed states do not accidentally have a tiny mass difference, 
one should be able to neglect the mixing between hadro-bottomonium states. This means that the 
1 hadro-bottomonium states (possibly with a core hb of higher excitation 
cleanly either into T(nS')7r7r or into hf)(mP)TT7r. 


23l|l. should decay 


III. TETRAQUARK 
A. Assumptions 


Tetraquarks are four-quark states constructed in analogy to the regular quark model. In par¬ 
ticular, the quarks are held together by effective gluon exchanges. Thus a necessary feature of 
each tetraquark model is that each isoscalar state is accompanied by (nearly) degenerate isovector 
states analogous to the p-u degeneracy in the light meson sector. 

Various variants for tetraquark models can be found in the literature. As one representative of 
this class of models we here discuss in detail only the implications of the most recently proposed 
interaction by Maiani et al. 3^ (see also the review article Ref. 37111. Note that the interaction 


originally proposed for the hidden charm states by the same group 


38l | was shown to be inconsis¬ 


tent with the most recent discussions [3a, [39(] . In this model tetraquarks are understood as such 


^ For both transitions, the dipion system is dominantly S-wave, and the available phase space for either of them is 
larger than 560 MeV for the tttt invariant mass. In this range, one expects that the /o(500), thus the tttt final state 
interaction (FSI), plays an important role. The FSI was neglected in Ref. [^. However, it was found in Ref. 
that the tttt FSI reduces the values of chromo-polarizability for the cc and bb by a similar factor. One should also 
keep in mind that the charged bottomonium-like Zj^(10610) and Z^{10650) states couple to both the Ttt^ and 
hence contribute to the dipion transition between T' and T. Such a contribution has not been taken into 
account in the analysis of the chromo-polarizability so far. 
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compact diquark-anti-diquark bound systems that the spin-spin interactions within the tetraquark 
are dominated by those within the diquarks. 

In this model, the mass of a tetraquark is given by 36|] 


M — Mqq - 2aL ■ S 2Kcq [(sg • Sc) -|- {sq ■ Sg)] , 


(9) 


where Sf{f = q,c,q,c) are the spins of (anti-)quarks, S is the total spin, L is the orbital angular 
momentum between the diquark and anti-diquark. The (anti-)quarks within the (anti-)diquarks are 
assumed to be in an S'-wave. The parameters Mqo, Be, a and Kcq are to be fixed from experiment. 
Denoting the spin of the diquark and anti-diquark as s = Sg -|- Sc and s = Sg -|- Sg, respectively, the 
Hamiltonian can be evaluated for a given tetraquark state of total angular momentum J, denoted 
by |s,s; S,L)j, 


M — Mqq Be —-— -- -|- Cl[T(T 1) S[S 1) — J{J 1)] -|- Keq [^(s -|- 1) -|- s{s -|- 1) — 3] . (10) 


For J = 1 the expression agree to Eq. (38) in Ref. 


36l |. Note that the parameters Be, a and k, 


are positive values, extracted from the experimental data by Maiani et al. 


cq 


36l |. As a result, the 


mass of the tetraquarks increases with increasing L and S, but decreases for growing J, which is 
a rather unusual feature. 


B. Consequences 


A general feature of tetraquark models is that a very rich spectroscopy emerges. In addition, 
there are always approximately degenerate isospin singlet and isospin triplet states, analogous to 
the case of the p and oj for the traditional qq mesons. 


Following Ref. 


36l | we will discuss the implications of the above model for S'-wave and P-wave 


tetraquark states only. The identiheation of some of the tetraquark levels with observed states was 
presented already in Ref. [36(. Here we extend this investigation by discussing all possible states 


with the mentioned quantum numbers. 

For S-wave tetraquarks, since L = 0 and J = S we can use |s,s)j to abbreviate \s,s; S,L) j. 
Then it follows from Eq. (IlOp that there are three sets of tetraquark states, whose masses are 
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4.2 

- 

(a) 





4.6 

(b) ^B630) 

4.0 




Z,(4020) 


> 

4.4 

— 

- 





_ _ 



X(3915) 

X(3872) 

- 9 - 

Z,C3900) 

—— 

X(3940) 

o 

!/5 


7(4290) 

- 

3.8 






C<3 

s 

4.2 

y(4260) 









r(4008) _ 

3.6 







4.0 




0++ 

T+ 

F- 

2 ** 



Q— p- 2" 3" 0“+ 2 -* ■ 


FIG. 2: The S'-wave (a) and P-wave (b) tetraquark spectroscopy in the charmonium sector. The green 
dashed and blue solid lines are input and prediction, respectively. 


approximately degenerate within the same set, since s and s are equal:® 

'o++: |0,0)o; 

< 1++: ^(|l,0)i + |0,l)i), 1+-: ^(|l,0)i-|0,l)i); (11) 

J++: ll,l)o, |1,1)2, 1+": 


For each of the quantum numbers, as mentioned before, the model predicts an isospin singlet 
as well as an isospin triplet.® Therefore, there should be in total 24 S'-wave tetraquark states 
already even without considering radial excitations that are expected about 400 MeV heavier than 
the corresponding ground states [36|. Among the ground states, the authors of Ref. [3^ identified 


the X(3872) as one of the neutral l"*"^ states and the Zc{3900) and Zc(4020) as the iso-triplet I"*" 
states, (|l,0)i — |0, l)i) /\/2 and |Ll)i, respectively. One may also assign A"(3915) and A(3940) 


as ll,l)o and | 1 , 1 ) 2 ) respectively 


36l |. although for each of them there is quite a large deviation 


between the mass of the tetraquark predicted and the actual mass of the observed state, cf. Fig. [2l 
Therefore, at least 15 of the 24 S'-wave tetraquarks are waiting for an observation. 

For P-wave tetraquarks, four isospin singlet 1 states without radial excitation were discussed 
in Ref. 3^ . But there are many more states — a few with exotic quantum numbers (that can not 


® The quantum numbers are , where the C-parity is given for the iso-singlet and the neutral member of the 
iso-triplet. 

® If we take into account the strange and anti-strange quarks, there would be a nonet for each . 














10 


be reached by the conventional qq states) like 0 and 1 

f 

1 " 

0 " 


|0,0;0,l)i; 
1 




(|l,0;l,l)o + |0,l;l,l)o), 


1 


^(| 1 , 0 ; 1,1)2 + 10 , 1 ; 1,1)2), 
|l,l;0,l)i; 

|i,i; 1,1)0; 

|l,l;l,l)i, 2”: |1,1;2,1)2; 

11 , 1 ; 1,1)2; 

|l,l;2,l)i; 

|i,1;2,i)3. 


V2 


(|1,0; 1,1)0-|0,1; 1,1)0); 


^(|1,0; 1,1)1-|0,1; 1,1)1); 


V2 


(| 1 , 0 ; 1,1)2-10,1; 1,1)2); 


( 12 ) 


According to Eq. (fTH]) the states in each line are (approximately) degenerate. In addition, as before, 
all of them appear as an isospin singlet and an isospin triplet when we restrict the light quark and 
anti-quark to be the up and down flavors. Thus, the above list amounts to 56 P-wave tetraquark 
states even without radial excitations taken into account. 

The assignments of the 1 P-wave tetraquark states to the observed structures in Ref. 36|] are 


as follows: among the four 1 states, three of them were identified with the Y (4008), Y (4260) and 
y(4630),^ and the other one was identified as one of two structures, called +(4220)® and +(4290) 


observed in e“’'e —^ hr'K'^TT 


4^. The states + (4360) and + (4660) were suggested to be the radial 


excitations of the +(4008) and +(4260). Thus, up to the first radial excitation, only 6 of 112 (28 
if considering only the isospin singlet ones) P-wave tetraquark states have candidates so far. 

There is one salient feature of Eq. m- for states with the same s, s, S and L, the mass 
decreases for increasing J, which appears to be a consequence of the negative sign in front of 
J( J + 1). Eor instance, among the states |1,1; 2,1) j with J , the 1 state has the largest mass 
while the 3 one has the smallest. Thus the observation of a rather light charmonium with J = 3 


could provide strong support for the tetraquark picture of Ref. 


361]. 


Besides the model discussed in detail above, also other tetraquark models can be found in the 


^ It was proposed in the literature that +(4630) observed in the AcAc and +(4660) observed in correspond 

to the same state [ 4 ^. 1^. 

® It was proposed in the literature that +(4220) observed hc'K'K and +(4260) observed in Jj'tp'K'K correspond to the 
same state [13, 
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literature that differ in the underlying assumptions. For example, in Refs. 4J-l48| the states are 
treated as four-quark systems without any clustering into diquark-anti-diquark assumed. As a 
result, the color part of the wave functions includes both antitriplet and triplet or sextet and anti¬ 
sextet conhgurations for the quark and antiquark pairs, respectively 


48l ]. As a result the number 


of the predicted S'-wave tetraquarks (Fig. 2 of Ref. [481]) is twice as large as that of Maiani et al. 
Although this picture can explain certain phenomena such as the narrow width of X (3872) due to 


its tiny J/il) + p and J/if: + uj component in the wave function [4^, Y (4140) as the hidden strange 
analog ccss of A(3872) 48], and Z(,(10610) and Z(,(10650) as bbqq four-quark systems Q, there 
is an even larger number of tetraquarks waiting to be observed within this approach. 


IV. HADRONIC MOLECULES 


A hadronic molecule is an extended object that results from non-perturbative scatterings of two 
or more hadrons. The hadronic molecules of interest here are bound states of a pair of charmed or 


bottomed mesons, which are similar to the deuteron as a bound state of the proton and neutron 


Q- 


Since the masses of some of the X, Y, Z states are close to S'-wave thresholds and couple strongly to 



Naively, one might expect that the number of possible molecules is at least as large as that 
of the available S-wave thresholds. In addition, since the open-charm and open-bottom mesons 
carry isospin 1/2 a pair of them can couple to both isospin 0 and 1. One might therefore expect 
almost degenerate isoscalar and isovector states for each quantum number similar to the tetraquark 
scenario. However, both expectations are not correct. First of all, a shallow bound state with an 
unstable constituent can in general not be narrower than that constituent, but will typically be 
broader 




59l |. In addition, the life time of a broad hadron, whose width is of the order of or 


even larger than the inverse of the range of forces, is too short to form a bound state with another 


hadron 


60l |. Thus, only the narrow L)i(2420) with a width of ~ 25 MeV can form an observable 


hadronic molecule (examples will be discussed below), while the broad Hi(2430) with a width of 
~ 380 MeV cannot. In this sense it also appears natural that the widths of y(4260) and V(4360) 


are of order 100 MeV. 
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In addition, the scattering potential in general comprises two contributions: a long-ranged part 
mediated by one-pion exchange and a short-ranged part that is often parametrized as contact 
interactions (and that one might be phenomenologically identihed with the exchange of heavier 
mesons, and could also come from s-channel qq states or more complicated dynamics). The short- 


ranged part needs to be fixed from data, as done in Refs. 


61I46,'1I| for the systems of a pair of 


S-wave heavy and anti-heavy mesons using the information of the X(3872) and .^^(lOGlO) states 
as input. Since at present there is not enough experimental information available for doing this 
for all channels, in this Section we have to restrict ourselves to qualitative statements regarding 
the hadronic molecular picture.® However, this has already allowed us to highlight some striking 
differences in the features of this model in comparison with the tetraquarks and hadro-charmonia. 

As a guidance for the existence of shallow bound hadronic molecules (here, we wish to use the 
phrase “shallow” in a loose sense meaning states with binding energies significantly less than 100 
MeV), we will study the contribution of the one-pion exchange and argue that if the one-pion 
exchange is repulsive for a given system the appearance of a bound state is unlikely, while a bound 
state could exist for an attractive one-pion exchange.^® This kind of argument is justified by the 
very small mass of the pion together with the observation that an attractive interaction mediated 
by a massless exchange particle always binds, regardless how weak the interaction is 
also notice the argument by Eriscon and Karl 


87|. We 


suggesting that two hadrons with an attractive 


one-pion exchange potential should form hadronic molecules if their reduced mass is sufficiently 
large. Based on this kind of reasoning it was possible to predict the existence of the X(3872) well 
before its observation There are also examples discussed below where the one-pion exchange 
does not contribute and possible molecular states must then be bound either by coupled channel 
effects or by short-ranged interactions. In such a case, the reasoning used in this paper cannot be 
applied. 

As mentioned above when talking about bound systems of mesons only those meson pairs are 
relevant where both mesons are sufficiently narrow. Here we focus on the ground state open charm 
meson doublet D and D* (characterized by a charm quark and a light anti-quark contribution 


® For the discussion of the one-pion exchange in effective field theories for the X(3872), see Refs. (Hll. Id^ - ldil ]. In 
addition, there exist various model calculations based on one-meson exchange potential or SU(4) extension of the 
light meson interactions. For the systems of a pair of S'-wave heavy mesons, see, e.g., Refs. |7CH83| . For the systems 
of an S-wave and a P-wave heavy meson, see Refs. [H^. . 

Note that in Ref. it is claimed that one-pion exchange does not contribute to the binding of, e.g., Zc(3900), 
since it gets cancelled by the contribution of the rj and q' exchanges in the P(3) limit. However, especially in the 
D* D system where the exchanged pion is near on-shell while q and q' are far off-shell, one should expect sizable 
violations of 17(3) symmetry. 
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(a) (6) (c) 

1/2 ffni) 1/2 1/2 ff^i) 1/2 1/2 5'2T(i) 3/2 

- 1 - - 1 - - 1 - 

I I I 

TT TT TT 

I I I 

_I_ I_ I_ 

1/2 ff T( 2 ) 1/2 3/2 ffi T( 2 ) 3/2 3/2 ff2 T( 2 ) 1/2 

FIG. 3: The various one-pion exchange diagrams relevant for the hadronic molecular scenario. Diagrams 
(a) and be show the t-channel pion exchange for the 1/2 + 1/2 and the 1/2 + 3/2 system, respectively, while 
diagram (c) shows the w-channel exchange, relevant for the 1/2 + 3/2 system only. 


with sf = where S£ is the total angular momentum of the light part) and the spin doublet 
that contains Hi (2420) and H2(2460) (characterized by a charm quark and a light anti-quark 
contribution with total sf = I"*") and their anti-particles as possible constituents. In particular, 
we discuss the ^ ^ and | +1 hadronic molecules, where we have used the total angular momentum 

of the light quark contribution to characterize the states. 

The Feynman diagram of the one-pion exchange is shown in Fig. [3l Because the pions are 
in the adjoint representation of the isospin SU{2) group, and the non-strange heavy mesons are 
isospin-1/2 states, each heavy-meson-pion vertex is accompanied by a factor r which stands for 
the Pauli matrices operating in the isospin space. The full one-pion exchange contribution between 
two heavy mesons therefore comes with a factor 7 / 1 ) • 7 / 2 ), where the subindices label the heavy 
meson to which the corresponding vertex is attached. One finds for a given pair of isospin-1/2 
particles with total isospin I and the third component 


■T(^2)\Ih) - 2 


^(^ + 1) 



(13) 


which is —3 for 1 = 0 and 1 for 1 = 1. Thus, if for a given set of quantum numbers the one-pion 
exchange contribution is attractive in the isoscalar channel, it will be repulsive in the isovector, 
and vise versa. Consequently, there is either an isoscalar state or an isovector state with some 
fixed quantum numbers, but typically not both. This is in contrast to the predictions for the 
tetraquark models discussed in the previous section, since the one-gluon exchange potential is 
flavor independent. 

In this paper we discuss where to expect S'-wave molecular states from non-perturbative in¬ 
teractions of the ^ multiplet {D,D*) with the anti-particles {D,D*) or with the I"*" multiplet 
(Hi, D 2 ). Particles and anti-particles need to be combined such that the states have a well defined 
C-parity. Bound systems formed from members of two | multiplets are potentially too broad to 
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show a striking signal in experiment and will therefore not be discussed here. Within this model 
space the following quantum numbers can be reached for ^ 


0 ++ : 
1 +- : 

1 ++ : 
2 ++ : 


\D*D*), \DD); 

-^\DD* - DD*), \D*D*)-, 

V2 

-^\DD* + DD*) 

v2 

\D*D*). 


(14) 


Here and in the following, we use the phase convention for the charge conjugation so that 
Dj = CDjC~^ with Dj representing any charmed meson and C the charge conjugation opera¬ 
tor. Whenever there appears more than one state for given quantum numbers, hadronic molecular 
states may appear as a result of coupled-channel dynamics and very limited statements are possi¬ 
ble without a detailed dynamical calculation. Exception to this are the I"’""*' and states. It is 
natural to identify the I"’""'' state with the X(3872). In addition, it turned out that HQSS forces 
the binding potentials for these two cases to be equal at LO in the low-energy expansion such that 
one arrives at the prediction of a tensor state X 2 located close to the D*D* threshold [63(]. In 
Ref. 90] it was shown that within a scheme where the one-pion exchange is treated perturbatively, 


this result is stable under inclusion of the DD inelastic channel in the D wave. As mentioned 
above, as a result of the isospin factor of the one-pion exchange, in the molecular picture one does 
not expect any isovector states with the quantum numbers 1^'*' and 2"’'+. 

Coupled-channel equations based on LO S-wave interactions typically produce as many states 
as the channels included in the calculation. Thus, we expect two states to be present with quantum 
numbers T' that strongly couple to DD* and D*D*. Indeed, both are established experimentally 
with the isovector states Zc(3900) and Zc(4020). Current experimental evidence locates both states 
above the thresholds For S-wave states this can be achieved either by momentum dependent 
interactions 9]| or non-trivial coupled-channel dynamics. Which one of these mechanisms is at 


work here (if any) requires a detailed model-building which however is beyond the scope of this 
work. Again, since isovector states are observed one should not expect their isoscalar partner states 
with 1^ quantum numbers. 

The situation for is more complicated, since also for this system we are faced with a coupled- 
channel system. In addition, the one-pion exchange is not allowed for the diagonal DD interaction 
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It should be stressed that the location of the poles is not settled yet: that the peaks in the data are located above 
the threshold does not necessarily imply that the poles are above the threshold as well. This is demonstrated for 
the ^6-states in Ref. [57| . 
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as a consequence of parity conservation. Neither does HQSS equalize the LO contact interaction in 
this channel to that in any other channel j^. Thus, at this point we cannot make any statement 
about the existence or non-existence of hadronic molecular states in the channel. 

The number of available channels and quantum numbers for the ^ + | system is even much 

2 2 " 

1 


larger than the one for h + h- 


1-+ : 

^|PPi+PiP), 

1,- 

1” : 

— iPPi-PiP), 

2"+ : 

■^\DD2+D2D), 

1,- 

2" : 

—\DD 2 -D 2 D), 

0"+ : 

-^\D*Di + DiD*) 

72 

0” : 

- DiD*) 

72 

3"+ : 

-^\D*D2 + D2D*) 

72 

1 , - 

3” : 

— |P*P2 - D 2 D*) 

72 


1 

V2 

1 

71 

1 

71 


\D*Di - DiD* 


D*Di - DiD*) 


-^\D*D2+D2D* 

V2 
1 


72 


D*D 2 - D2D* 


^\D*D2+D2D* 

v2 


72 


D*D2 - D2D* 


(15) 


In this case, one can also check in which channels HQSS predicts the same LO interaction. For 
each isospin, 0 or 1, there are four independent interactions denoted as (s^i, S£2) 72 ’ 

+ ^7- 


with sl = sei + se 2 = sL + sL: 


^ ^ 1 
2 ’ 2 ’ 


Vr 


^ ^ 1 
2 ’ 2 ’ 


^ ^ 1 
2 ’ 2 ’ 


Vr 


^ ^ 1 
2 ’ 2 ’ 


1.1,2 

2 ’ 2 ’ 


Vr 


i ^ 2 
2’ 2’^ 


i ^ 2 
2’ 2’^ 


Vr 


^’^’2).(16) 


It turns out that for the diagonal interactions all of the above listed channels have a different linear 
combination of the matrix elements given in Eq. (jl6p . Therefore, not much can be derived from 
HQSS without further input in this case. 

In addition, in case of the ^ + | system two kinds of one-pion exchange contributions are 
possible in the channels with odd parity. They correspond to the t-channel and u-channel one- 
pion exchange, cf. Fig. [3l Denoting the coupling constants for the coupling of a pion with two 
I states (P-wave), with two | states (P-wave) and connecting a ^ state and a | state (P-wave) 
by g, gi and 52 , respectively, the potential for the t-channel exchange is proportional to ggi and 
the one for the u-channel exchange is proportional to g^- While the magnitude and sign of the 
coupling constant g are fixed by data and lattice QCD calculations, nothing is known yet about gi. 
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We thank J. Nieves for pointing this out. 
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FIG. 4: The two-body thresholds in the charmonium mass range potentially relevant for the formation of 
hadronic molecules. 


Therefore, we cannot identify the channels where the corresponding one-pion exchange potential is 
attractive. The situation is different for the part of u-channel exchange with two D-wave vertices 
proportional to (although the additional momentum factors on the vertices might suppress this 
part of the potential). Looking only at the sign of this tt-channel exchange contribution for the 
uncoupled channels, one finds attraction for isoscalar states with the exotic quantum numbers 0 
and 3 while this part of the interaction is attractive in the isovector channel with the quantum 
numbers 0 and 3 . However, as mentioned above, since for these systems we cannot even 

analyze the complete one-pion exchange contribution no strong conclusion can be drawn regarding 
the existence or non-existence of these states. 

For the quantum numbers 1~^, 2“^ three channels couple to each other. Thus, here even less 
can be said without dynamical analysis. However, some general statements are still possible. For 
instance, since hadronic molecular states appear in the vicinity of S-wave thresholds, the lightest 
molecular state in the 5 + | family is expected to be close to the lowest threshold, namely__DiL> 
(cf. Fig. 01), and with the quantum numbers of either 1 or 1 Indeed, in Refs. |52l. l53l| the 


M 


Y (4260) was proposed to be (predominantly) a DiD bound state with = I , and in Ref. 
the molecular states with exotic quantum numbers = I were investigated. 

For the quantum numbers 0“^, if a corresponding state exists at all, it should be around the 
D*Di threshold and thus more than 100 MeV above the mass of the y(4260). While this mass 
range is similar to what is predicted in the tetraquark picture, it is very different to what is 
predicted in the hadro-charmonium scenario, where the 0 states are predicted to be the lightest 
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of their kind. 

One striking distinction between the tetraquark picture and the hadronic molecular scenario 
becomes visible for the J = 3 states; while in the tetraquark model of Ref. 3^ those are amongst 


the lightest states, in the hadronic molecular picture, if they exist, they should be close to the 
D 2 D* threshold, again more than 100 MeV above the y(4260), cf. Fig. [H 

V. SUMMARY 


In this work we have investigated the exotic charmonium spectrum in different scenarios based 
on HQSS. In particular, we have compared the spectra arising from three models, i.e. hadro- 
charmonium, tetraquark and hadronic molecule, which turn out to provide quite distinct predictions 
for not yet discovered “AYZ” states. 

In the hadro-charmonium scenario we find that if y(4260) and y(4360) are mixed hadro- 
charmonium states with = 1 as proposed in Ref. 27] then two spin partners with = 
0 and two more with = 1 ^ and = 2 respectively, should exist as well. Their 
possible production and decay modes are discussed to guide the experimental search in the future. 

In the tetraquark scenario the spectrum is much richer than in other scenarios with approxi¬ 
mate degeneracies between isospin singlet and triplet states. If one assigns the observed A(3872), 
Zc(3900) and Zc(4020) to the 5-wave states, some of the states with = 1 can also be 
assigned to the P-wecve tetraquark states, i.e. y(4008), y(4260) and y(4630), while the possible 
existence of two 1 states, y(4220) and y(4290) in e"*“e“ —)• hcTT~^TT~ will accommodate all the P- 
wave 1 tetraquark ground states. However, it should be noted that so far only a small number of 
predicted tetraquark states can be assigned to existing observations — at least 67 among 80 of the 
ground states are left to be discovered, if compact tetraquarks provide the dominant components 
of the new XYZ states. 

In the scenario of hadronic molecule very little can be said without detailed modeling of the 
relevant hadron-hadron potentials. Under the assumption that the one-pion exchange potential 
plays a dominant role, one expects bound states to appear either in the isoscalar or isovector 
channel but not in both simultaneously. The observation of the isospin singlet A(3872) and the 
isospin triplet Zc(3900) and Zc(4020) without evidence of isospin partners matches this expectation 
well. Taking into account that the constituent mesons are to be narrow enough in the hadronic 
molecular system, the allowed number for hadronic molecules formed by low-lying narrow charmed 
meson pairs turns out to be a lot smaller than what is predicted within the tetraquark scenario. 
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In order to highlight that the spectroscopy of the different scenarios is indeed very different 
we summarize the implications for two quantum numbers: Amongst the three models only the 
hadro-charmonium scenario predicts 2 0 states with one of them being even lighter than the 
F(4260). In the tetraquark picture also two 0 ^ states appear, however, both are predicted to 
be heavier than y(4260). Due to lack of a dynamical model, in the hadronic molecular picture 
there is no prediction yet for a 0 '' state, however, if it exists it should be located near the DiD* 
threshold — again above the mass of the y(4260). 


Interestingly, we find that the tetraquark model of Ref. 


36l | and hadronic molecular scenarios 


predict strikingly different patterns for the J = 3 states. In the former scenario the = 3 
state is expected to be among the lightest states, while in the hadronic molecular scenario the 
J = 3 states should have masses close to the D 2 D* threshold if they exist. An observation of the 
J = 3 states should provide crucial information about the underlying dynamics. 

It is very important to test the different scenarios in future experimental studies at LHCb, 
BESIII, Belle-II, PANDA and others. 
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